Quantum Necking in Stressed Metallic Nanowires 
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When a macroscopic metallic wire is subject to tensile stress, it necks down smoothly as it elon- 
gates. We show that nanowires with radii comparable to the Fermi wavelength display remarkably 
different behavior. Using concepts from fluid dynamics, a PDE for nanowire shape evolution is 
derived from a semiclassical energy functional that includes electron-shell effects. A rich dynamics 
involving movement and interaction of kinks connecting locally stable radii is found, and a new class 
of universal equilibrium shapes is predicted. 
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It has recently become possible to image metal 
nanowires with subatomic resolution, and record their 
structural evolution in real time, using transmission- 
electron microscopy 0, 0, 0, S 01 1 thus opening the door 
to a world of nonlinear dynamical phenomena at the 
nanoscale. Heretofore, most information about the struc- 
ture and dynamics of nanowires was obtained indirectly, 
from transport and cohesive properties la M. la , |9l llCfl . 
Classical molecular dynamics simulations [ll|, |TJ, |ljj, |b| 
provide important insights into atomic structure and dy- 
namics of nanowires, but cannot explain the observed 
stability Q, 0, |j| of long nanowires, which would be un- 
stable classically under surface tension [l5| . Quantum 
molecular dynamics can in principle study this stability, 
but has so far been limited to systems too small to ad- 
dress such issues (e.g. single or double chains (lH.ll7| ). 

By contrast, a nanoscale free-electron model 0, 0, 
I20L Elj has proven successful in explaining the quantum 
suppression of the Rayleigh instability [20(. Here we ex- 
tend this approach by developing a dynamics for a free- 
electron nanowire under the assumption that surface dif- 
fusion is a dominant process. This dynamics, of instrinsic 
interest in nonlinear science, explicitly includes quantum- 
size effects and allows the study of a broad range of dy- 
namical processes (see Figs. ^|3J|: Approach to equilib- 
rium of a nanowire, propagation of an instability, and 
evolution of a wire under elongation/compression. Fig- 
ure^shows how an initially random nanowire necks down 
to a universal shape, consisting of a central cylinder con- 
nected to contacts having the form of unduloids of revolu- 
tion, providing a dynamical mechanism for the nanofab- 
rication technique invented by Kondo and Takayanagi 
0. Figure compares well to early experiments on gold 
nanocontacts showing perfect correlation between 

force and conductance with steep conductance steps, and 
a hysteresis between elongation and compression. 

The model of a nanowire consists of free electrons con- 
fined to an axisymmetric wire by a hard-wall potential 
[l8| . The wire radius is R(z,t) in the interval [0 , L] 
along the z-axis, and periodic boundary conditions |22( 
are used to extend R(z). In contrast to the classical 



case for which all nonaxisymmetric distortions raise 
the energy, Jahn- Teller distortions that break axisym- 
metry of nanowires do lower the energy in a few cases. 
However, the deepest energetic minima occur for axisym- 
metric wires |23| - and axisymmetry, once present, is pre- 
served by the diffusion equation, so that Jahn- Teller dis- 
tortions are most likely suppressed dynamically. 

In the spirit of the Born-Oppenheimer approximation, 
the total energy of the nanowire is taken to be the elec- 
tronic energy. Since we are dealing with an open system 
of electrons, the grand-canonical potential is used, and 
can be separated into Weyl and mesoscopic contributions 
|l9l I24L Eq . Dropping the volume contribution, which is 





FIG. 1: (a) Electron-shell potential V(T, R) for T = 0.002T F . 
Here ef and Xf are the Fermi energy and wavelength, Tf = 
EF/kB and kF = 2tz/Xf- (b)-(d) Approach to equilibrium of 
an initially random wire: (b) initial shape; (c) r = 2 x 10 4 ; 
(d) t — 3 x 10 7 , equilibrium wire with G = 12 Go. 



2 



assumed to be constant, the energy functional is 

a[T,R(z)]=a(T)S[R(z)]+ [ V(T,R(z))dz, (1) 

Jo 

where cr(T) is the surface tension [2||, S is the surface 
area of the wire, and V is a mesoscopic electron-shell 
potential whose integral gives the so-called "shell correc- 
tion" to the energy (as used, for instance, in the study of 
clusters H3)- Higher-order terms 0, proportional 
to the mean curvature, etc., are neglected. V can be 
expressed in terms of a Gutzwiller-type trace formula 

V(T,R) = ^ E E > ( 2 ) 

W—l V—ZW 

where the sum includes all classical periodic orbits (v, w) 
in a disk billiard j^, characterized by their number of 
vertices v and winding number w, L vw = 2vRsui(ttw/v) 
is the length of orbit (v,w), and 9 VW — kpL vw — 3vir/2. 
The factor f vw — 1 for v — 2w, 2 otherwise, accounts for 
the invariance under time-reversal symmetry of some or- 
bits, and a vw (T) = t vw / sinhr^ (t vw = nk F L vw T/2T F ) 
is a temperature-dependent damping factor. A similar 
trace formula was previously derived [20j for small ax- 
isymmetric perturbations of a c ylin der using semiclassi- 
cal perturbation theory (23, |23, |2jJ . Here we point out 
that semiclassical perturbation theory remains valid for 
large deformations, as long as new classes of non-planar 
orbits can be neglected (adiabatic approximation). 

Eqs. and J5J define an energy functional which 
is general and simple enough to solve for nontrivial 
nanowire geometries with a wide range of radii. The mul- 
tiple deep minima of V(R), shown in Fig. Ufa), favor cer- 
tain magic radii and suggest a multistable field theory 
analogous to the sine-Gordon model. We therefore an- 
ticipate solutions consisting of cylindrical segments con- 
nected by kinks. 

Rather than directly minimizing the energy func- 
tional |Jy , we derive an equation describing the nanowire 
surface dynamics which yields the approach to equilib- 
rium as well as the onset of instability. The model has 
two main assumptions. First, in the thin wires con- 
sidered, the majority of atoms are at the surface and 
thus surface self-diffusion is the mechanism of ionic mo- 
tion. Second, under the Born-Oppenheimer approxima- 
tion, the electronic energy acts as a potential for the 
ions, and the dynamics derive from ionic mass conserva- 
tion: 

8R{z,t) v a d 



dt 



R(z,t)dz^ J ^ (3) 

where v a = iir 2 /k F is the volume of an atom, and the 
z-component of the surface current is given by Fick's law: 
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where R 1 = dR/dz, and ps and D$ are the surface 
density of ions and the surface self-diffusion coefficient, 
respectively. The precise value of Dg for alkali met- 
als is not known, but it can be removed from the evo- 
lution equation by rescaling time to the dimcnsionlcss 
variable t = (psDsTp/T)t. The chemical potential 
fi is obtained by calculating the change in the energy 
(|T|) with the addition of an atom at point zq, p(zq) = 
Cl[T,R(z) + AS(z-z )] -n[T,R(z)], where A= v a j2-KR 
is chosen so that the volume of an atom is added: 
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where b vw (T) = a vw cosh r vw . Eq. © is then solved nu- 
merically using an implicit scheme. Our non-linear dy- 
namical model, Eqs. J2HHJ), differs from previous studies 
of axisymmetric surface self-diffusion by the 

inclusion of electron-shell effects [2nd line of Eq. (JSJ)], 
which fundamentally alter the dynamics. 

To explore the equilibrium shapes of nanowires, we 
considered several random initial configurations of var- 
ious lengths and widths, one of which is shown in Fig. 
^b). Under surface self-diffusion, the short- wavelength 
fluctuations are rapidly smoothed out [Fig. [He)], an d 
several cylindrical segments connected by "kinks" are 
visible. On a much longer time-scale, the wires neck 
down to a universal equilibrium shape consisting of a 
nearly perfect cylindrical region with a radius near one 
of the minima of V(R), connected to thicker leads [Fig. 
IHd)]. In all fifteen cases studied, the leads closely ap- 
proximate Delaunay unduloids of revolution. The undu- 
loid is a stationary configuration, unstable under surface 
tension alone |22j, but stabilized in the present model 
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FIG. 2: Propagation of an instability for a wire of initial 
radius /cf-R = 10, time progressing from bottom to top (At = 
15), with curves shifted for clarity. The top curve corresponds 
to a much later stage in the time evolution (r = 6 x 10 4 ). 
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FIG. 3: Shape of a nanowire connecting two quasi-macroscopic leads during elongation (left column, a-g) and subsequent 
compression (right column, h-n) p^l . For comparison, the entire system of length L = 280 kp 1 + AL is shown in inset (i), 
the small square indicating the nanowire shown in image (a). Adjacent wires have the same length, from top to bottom: 
kpAL — 0, 10, 20, 40, 70, 105, 140. Above each wire is its conductance in units of Go. Central Inset: Conductance and force as 
a function of AL: (i) Conductance during elongation (solid line) and compression (dotted line); (ii) force during elongation; 
(iii) force during compression. 



by electron-shell effects in the connecting cylinder. The 
unduloids arise as superpositions of large numbers of 
closely spaced kinks. Equilibrium wires with conduc- 
tance G/Gq = 1, 3, 6, 12, 17, 23, . . . were obtained, where 
Go = 2e 2 /h is the conductance quantum. Intermediate 
conductance values may occur more rarely. 

Other cases of interest are the propagation of an insta- 
bility [3^ | and the homogeneous-inhomogeneous transi- 
tion in a nanowire [2l| . In the former case, starting with 
a cylinder of an unstable radius, a maximum in V(R) 
[Fig.^a)], and adding a localized Gaussian deformation 
of small amplitude (0.001.R), we find that the perturba- 
tion decreases during an incubation period and then the 
instability sets in with essentially a single Fourier com- 
ponent corresponding to the most unstable mode [2lj| . 
saturates at a finite amplitude, and propagates with con- 
stant velocity (see Fig. A similar instability is ob- 
served in the pearling of membranes pjjj ]. The classi- 
cal Rayleigh instability itself [34| propagates in much 
the same way, but does not saturate, instead leading to 
breakup of the cylinder. Both the wavelength [2lJ and 
front velocity of the instability are strongly dependent on 
the initial radius. Once the instability is fully developed, 
kink-antikink pairs start annihilating, eventually leading 



to a shape consisting of a series of cylinders of different 
radii [corresponding to neighboring minima in V(R), Fig. 
Cla)], connected by kinks. Starting from a more realistic 
initial condition with random perturbations, the instabil- 
ity occurs at several places more or less simultaneously. 

The dynamical model also allows for the study of the 
evolution of a nanowire under strain. The length L of the 
wire is changed at an average rate dL/dt in a sequence of 
infinitesimal steps. In each step, the wire deforms elas- 
tically, i.e. the change of length is distributed along the 
wire as if each of its slices were Hookean, with spring con- 
stant proportional to its cross-sectional area. The radius 
of each slice is simultaneously changed so as to conserve 
its volume. Between steps, the wire evolves according 
to Eq. 1(5}. We start with an equilibrated unduloid of 
outer radius fcpi? = 99, representing a nanocontact be- 
tween two quasi- macroscopic leads [see Fig. |3J (inset)]. 
We then strain the wire by elongating it at a constant 
speed kpdL/dr = 0.1 until it necks down to a long, 
atomically thin wire [Fig. OJb-g)], at which point we let 
it partially equilibrate [Fig.^h)]. Subsequently, we com- 
press the wire at the same rate until it reaches its initial 
length [Fig. Efi-n)]. During elongation, the necking of 
the wire occurs by nucleation of kink-antikink pairs at 
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the center with subsequent motion of these kinks toward 
the leads, leaving a long, cylindrical wire in the middle. 
During compression, the remaining kink-antikink pairs 
annihilate, but new kinks do not emerge from the leads. 
Similar necking was observed for various elongation rates 
kpdL/dr < 1. For much faster rates, structural relax- 
ation due to surface diffusion is suppressed. 

These different processes lead to a substantial hystere- 
sis between elongation and compression, as can be seen 
by computing the electrical conductance of the contact 
and the force applied to it. The force is obtained from the 
change in the energy Q with elongation, F = —dVl/dL. 
The conductance is computed from the Landauer for- 
mula, using the adiabatic and WKB approximations 0] 
to compute the transmission probabilities. Both quanti- 
ties are shown in the central inset of Fig. [31 Compared 
to our previous results obtained without structural relax- 
ation |18|. the conductance steps and relaxations of the 
force are much steeper, in better agreement with experi- 
ments 0, and a substantial hysteresis is observed in the 
conductance. It has been argued 0, Q that the abrupt- 
ness of the conductance steps, and their perfect correla- 
tion with the jumps in the force, rule out an electronic 
mechanism. However, we have shown that this behavior 
arises naturally in a model which takes proper account 
of electronic quantum-size effects on the structure, even 
when atomistic effects are neglected. 

A final interesting feature is the spikes in the force 
that occur at the opening and closing of a channel, cor- 
responding to the rapid energy relaxation during the 
creation or annihilation of a kink/antikink pair. These 
spikes are suppressed, and the hysteresis decreases, when 
the speed of deformation is increased, suggesting that the 
spikes might be observable experimentally by decreasing 
the elongation rate. A failure of this model is the force 
during compression: contrary to experiment 0, Q , the 
force is attractive, suggesting that a resistance to com- 
pression due to the ions needs to be added to the free- 
electron model. Indeed, an interesting question raised 
by this work is how discrete positive ions would arrange 
themselves to accomodate the predicted shapes favored 
by the conduction electrons. 

This research was supported in part by NSF grants 
DMR0072703 and DMR0312028 (CAS & JB), and by 
NSF grant DMR98 12526 (REG). JB also acknowledges 
support from the Swiss National Science Foundation and 
grant SFB276 of the Deutsche Forschungsgemeinschaft. 



[1] Y. Kondo and K. Takayanagi, Phys. Rev. Lett. 79, 3455 
(1997). 

[2] T. Kizuka, Phys. Rev. Lett. 81, 4448 (1998). 

[3] Y. Kondo and K. Takayanagi, Science 289, 606 (2000). 



[4 
[5 
[6 
[7 
[8 
[9 
[10 
[11 
[12 
[13 
[14 

[is; 

[16 

[17; 
[is: 

[19 
[20 

[21 

[22 

[23 

[24 
[25 

[26 

[27 

[28 

[29 
[30 

[31 
[32 

[33 
[34 

[35] 



V. Rodrigues and D. Ugarte, Phys. Rev. B 63, 073405 
(2001). 

V. Rodrigues and D. Ugarte, Phys. Stat. Sol. (b) 230, 
475 (2002). 

G. Rubio, N. Agrai't, and S. Vieira, Phys. Rev. Lett. 76, 
2302 (1996). 

A. Stalder and U. Diirig, Appl. Phys. Lett. 68, 637 
(1996). 

C. Untiedt, G. Rubio, S. Vieira, and N. Agrai't, Phys. 
Rev. B 56, 2154 (1997). 

A. I. Yanson, I. K. Yanson, and J. M. van Ruitenbeek, 
Nature 400, 144 (1999). 

A. I. Yanson, I. K. Yanson, and J. M. van Ruitenbeek, 
Phys. Rev. Lett. 84, 5832 (2000). 

U. Landman, W. D. Luedtke, N. A. Burnham, and R. J. 
Colton, Science 248, 454 (1990). 

T. N. Todorov and A. P. Sutton, Phys. Rev. B 54, 
R14234 (1996). 

M. R. S0rensen, M. Brandbyge, and K. W. Jacobsen, 
Phys. Rev. B 57, 3283 (1998). 

O. Giilseren, F. Ecrolessi, and E. Tosatti, Phys. Rev. 
Lett. 80, 3775 (1998). 

S. Chandrasekhar, Hydrodynamic and Hydromagnetic 
Stability (Dover, New York, 1981), pp. 515-74. 

H. Hakkinen, R. Barnett, A. Scherbakov, and U. Land- 
man, J. Phys. Chem. B 104, 9063 (2000). 

S. Bahn and K. Jacobsen, Phys. Rev. Lett. 87, 266101 
(2001). 

C. A. Stafford, D. Baeriswyl, and J. Biirki, Phys. Rev. 
Lett. 79, 2863 (1997). 

C. A. Stafford, F. Kassubek, J. Biirki, and H. Grabert, 
Phys. Rev. Lett. 83, 4836 (1999). 

F. Kassubek, C. A. Stafford, H. Grabert, and R. E. Gold- 
stein, Nonlinearity 14, 167 (2001). 

C. -H. Zhang, F. Kassubek, and C. A. Stafford, cond- 
mat/0209601. 

The locality of the semi-classical functional, Eq. sup- 
presses the dependence on boundary conditions. 

D. Urban, C.-H. Zhang, J. Biirki, H. Grabert, and C. A. 
Stafford, in preparation. 

V. Strutinsky, Nucl. Phys. A 122, 1 (1968). 
M. Brack and R. K. Bhaduri, Semiclassical Physics 
(Addison- Wesley, Reading, MA, 1997), and refs. therein. 
The value a — £Fkp/80ir was used |19l to model alkali 
metals. For noble metals, a is roughly twice as large. 

C. Yannouleas and U. Landman, Phys. Rev. B 51, 1902 
(1995). 

D. Ullmo, M. Grinberg, and S. Tomsovic, Phys. Rev. E 
54, 136 (1996). 

S. C. Creagh, Ann. Phys. 248, 60 (1996). 

B. D. Coleman, R. S. Falk, and M. Moakher, Phys. D 
89, 123 (1995). 

J. Eggers, Phys. Rev. Lett. 80, 2634 (1998). 

A. Bernoff, A. Bertozzi, and T. Witelski, J. Stat. Phys. 

93, 725 (1998). 

See physics . arizona. edu/~ st afford/necking. html. 

D. Zhang, T. Powers, R. Goldstein, and H. Stone, Phys. 
Fluids 10, 1052 (1998). 

R. Bar-Ziv and E. Moses, Phys. Rev. Lett. 73, 1392 
(1994). 



